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"cs" ! Abstract 

Within the framework of perturbative QCD approach, we study the charmless two-body 
decays B s — > ai(1260)(6i(1235))P(V) (P, V represent the light pseudo-scalar and vector mesons, 
respectively). Using the decays constants and the light-cone distribution amplitudes for these 
mesons derived from the QCD sum rule method, we find the following results: (a) The decays 
B® —?• a± K + (K* + ) have the contributions from the factorization emission diagrams with a large 
Wilson coefficient C2 + Ci/3 (order of 1), so they have the largest branching ratios and arrive 
at 1(T 5 order. While for the decays B° -> a?K°(K*°), the Wilson coefficient is C\ + C 2 /3 
in tree level and color suppressed, so their branching ratios are small and fall in the order of 
1CP 7 ~ 10~ 8 . For the decays B a — > b\K(K*), all of their branching ratios are of order few times 
1CT 6 . (b) For the pure annihilation type decays B® — > a\{b\)p except the decays B® —> aiir 
having large branching ratios of order few times 10 -6 , the most other decays have the branching 
ratios of 10 -7 order. The branching ratios of the decays B® — > o5(&5)w are the smallest and fall 
in the order of 10 -8 ~ 10~ 9 . (c)The branching ratios and the direct CP-asymmetries of decays 
^ , B® — > ai(bi)r]('' are very sensitive to take different Gegenbauer moments for ry- 1 '. (d) Except for 

the decays B® — > a\K*° , a^w, b\uj, the longitudinal polarization fractions of other B® — > ai(b\)V 

Q\ ' decays are very large and more than 90%. (e) Compared with decays B® — > a\{b\)P, most of 

W) • Bg — y aiib\)V decays have smaller direct CP asymmetries. 
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I. INTRODUCTION 

In general, the mesons are classified in J PC multiplets. There are two types of 
orbitally excited axial- vector mesons, namely 1 ++ and l + ~. The former includes 
ai(1260), /i(1285), /i(1420) and K\a, which compose the 3 Pi-nonet, and the latter in- 
cludes &i(1235),/ti(1170), /ii(1380) and K 1B , which compose the 1 P 1 -nonet. There is an 
important character for these axial- vector mesons except ax (1260) and b\ (1235), that is 
each different flavor state can mix with one another, which comes from the other nonet 
meson or the same nonet one. 

B° — > a^(1260)7r T are the first decay modes with an axial- vector in the final state 
observed by BarBar and Belle |l|-|3j. Measuring their time- dependent CP asymmetries 
can provide the information of Cabibbo-Kobayshi-Maskawa (CKM) weak phase a. After 
these measurements, many other charmless decays B — > AP, AV (P, V stand for the 
light pseudo-scalar and vector mesons) have also been reported by experiments |4J-|lO 
On the theoretical side, many met 
the naive factorization approach 
QCD factorization approach 



ichfn, 



rods are employed to research these decays, such as 



12] , the generalized factorization approach [13J, the 
, the PQCD approach [l6[ . Though the factorization 
approach holds only approximately and its predictions are at odds with experiments for 
some decays, many factorization approaches can explain the data in many cases. So 
these results predicted by the different factorization approaches are useful to investigate 
production mechanism of axial vectors in B meson decays, extract the information of 
Cabibbo-Kobayshi-Maskawa (CKM) weak phase, probe the structures of axial vectors, 
even calculate the relative strong phase between tree and penguin diagrams. To our 
knowledge there is still lacking the study of charmless decays B s — > AP, AV both in 
experiments and theories. Our aim is to fill in this gap and provide a ready reference to 
the forthcoming experiments to compare their data with the predictions in the PQCD 
approach. In view of the fact that a!(1260) and &i(1235) can not mix with each other 
because of the opposite C-parities and they do not also mix with other mesons, we would 
like to study the decays B s ->• ai(1260)P(V),&i(1235)P(V) in detail. 

In the following, oi(1260) and &i(1235) are denoted as a\ and b\ in some places for 
convenience. The layout of this paper is as follows. In SecJITl decay constants and light- 
cone distribution amplitudes of the relevant mesons are introduced. In SecJTTH we then 
analyze these decay channels using the PQCD approach. The numerical results and the 
discussions are given in Sec. IIVI The conclusions are presented in the final part. 

II. DECAY CONSTANTS AND DISTRIBUTION AMPLITUDES 

For the wave function of the heavy B s meson, we take 

®B 3 {x,b) = -==(p Bs +m B J 7 50 Bs (x,6). (1) 

Here only the contribution of Lorentz structure 4>b s {x, b) is taken into account, since the 
contribution of the second Lorentz structure <pB s is numerically small [17J and has been 
neglected. For the distribution amplitude cj) Bg (x,b) in Eq. ([!]), we adopt the following 



model: 



M 2 x 2 1 
fo.&b) = N Bs x 2 (l - x) 2 exp[-^y - -M) 2 ], 



(2) 



where Ub is a free parameter, we take Ub = 0.5 ± 0.05 Gev in numerical calculations, and 
Nb s = 63.671 is the normalization factor for ujb = 0.4. 

The wave functions for the pseudo-scalar (P) mesons K, n are given as 



^p(P,x,C) = 



V2N< 



=75 [lf(j) A (x) + m o P (x) + (mo (^ - v ■ n)<f) T (x)] 



(3) 



c 



where the parameter ( is either +1 or —1 depending on the assignment of the momentum 



fraction x. The chiral scale parameter m is defined as m 
amplitudes are expanded as: 

3fK,n 



Ml 



m qi +m q2 



- — . The distribution 



<j? K (x) 
<j) T K {x) 



3/k 

-Ik 

2^/Q 



3/2, 



,3/2, 



X(l - X) 1 + ai(K,7r)Cl it) + a2(K,7r)C 2 (*) 



l + 0.43C 2 1/2 (t) ;0P(x) 



2^ 



1 + 0.24C 2 1/2 (t) 



yl/2 



C7^(t) + o.35cr(0 ;CW 



1/2, 



-u 



2VQ 



yl/2 



C 1 i/ '(t) + 0.55C 3 i/z (t) 



1/2, 



with Gegenbauer polynomials defined as: 

Cl /2 {t) = 3t,C 3 2 /2 (t) = 1.5(5t 2 - 1), 

C\'\t) = t, Cl /2 (t) = 0.5(3t 2 - 1), ClJ 2 (t) = 0.5£(5£ 2 - 3). 



(4) 
(5) 
(6) 

(7) 



As for the distribution amplitudes of the pseudo-scalar mesons rj and rj', we use the quark 
flavor basis mixing mechanism proposed by Refs. 18( and take the same formulae and 
parameter values as those in Ref.|19|. 

For the vector mesons, their distribution amplitudes are defined as 



(V(P,el)\q 2 p(z)q la (0)\0) 
(V(P,e* T )\q 2 p(z)q la (0)\0) 



V2N C J 



dxe lxp - z [m v f L (l) V {x) + f L tf 
dx e lxp - z [m v 4* T 4> v v {x) + f T i 



(x) +m V (fi s v (x)} a p, 



V2N C J 
+m v ie llupa -f 5 ^e* T v n p v° '<t>y(x)] af) 



>v 



x 



(9) 



where n(v ) is the unit vector having the same (opposite) direction with the moving of the 
vector meson and x is the momentum fraction of q 2 quark. The distribution amplitudes 
of the axial-vectors have the same format as those of the vectors except the factor 275 
from the left hand: 



(A(P,e* L )\q 20 (z)q la (O)\O) 



(A(P,e* T )\q 2 p(z)q la (0)\0) 



«75 



V2N< 



dx e ixp - z [m A f L <j) A {x) + f L P<p A {x) + m A (f) A (x)} a p, 



c JO 



11$ 



dx e^ z [m A f T ^ A {x) + 

\/Zl\ c Jo 

+m A ie flupal5 Ye* T v n p v a (j) A (x)} a(3 . 



IX 



(10) 



TABLE I: Decay constants and Gegenbauer moments for each meson (in MeV). The values are 
taken at \i = 1 GeV. 



Ik* 

209 ±2 


Ik* 
165 ±9 


u 

231 ±4 


186 ±9 


Ik 
160 


130 


lai 

238 ± 10 


f T 

lb! 

-180 ±8 


209 ±2 


Jp 

165 ±9 


195 ±3 


f T 
151 ±9 


a\K 
0.17 





0.2 


0.44 


4m 

0.03 ± 0.02 


aim 

0.04 ± 0.03 


4m 

0.11 ±0.09 


4m 

0.10 ±0.08 


a|(p,w) 
0.15 ±0.07 


a£(p,u) 
0.14 ±0.06 


4{<i>) 

0.18 ±0.08 


4(0) 

0.14 ±0.07 


a|(ai(1260)) 
-0.02 ±0.02 


a] L (a 1 (1260)) 
-1.04 ±0.34 


oj[(6i(1235)) 

-1.95 ±0.35 


4(&i(1235)) 
0.03 ±0.19 



As for the upper twist- 2 and twist-3 distribution functions of the final state mesons, 
<fiv(A), 4>via)' < Pv(A)i ^v(AY 'fiv(A) anc ^ ^v(A) can ^ e calculated by using the light-cone QCD 
sum rule. We list the distribution functions of the vector (V) mesons, namely p(u, 0), as 
follows 
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The axial- vector (A) mesons , here a\ and bx, can be obtained by replacing each 0y with 
4>a, by replacing fy(fv) with / in Eq.ffTTT). Here we use / to present both longitudinally 
and transversely polarized mesons ax(bx) by assuming fj — f ai — f for ax and /^ = 
f b T = f for bx- In Eg. fill I) , the twist- 2 distribution functions are in the first line and can 
be expanded as 



6x(l — x) 



l + al 



2 V 



1) 



6x(l 



x) 



+ 3a}' t + a\' -(5t 2 -l) 



for V mesons; 



for A mesons, 



(12) 
(13) 



where the zeroth Gegenbauer moments a^a^ = %(bx) = and %(ax) = Q>o(bx 
As for twist-3 LCDAs, we use the asymptotic forms for V mesons: 



1. 



h 



(*). 



x 



3t 2 ,h 



W, 



x 



4l\x) 



6x(l — x) 



6x(l 
,gy{x) = 



X 



;i±t 2 ). 



(14) 



And we use the following forms for A mesons: 

h^(x) = 3a^t 2 + -aj;t(3t 2 - 1), h^\x) = 6x(l - x)(af + aft), 

g { l\x) = 6x(l-x)(al + alt),g<?\x) = hl(l + t 2 ) + h\t\ (15) 

In Eqs.tjl])-® and Eqs. ffT2|) - ffT5|) . the function t — 2x — 1. The decays constants and the 
Gegenbauer moments all' for each meson are quoted the numerical results 2CH25J and 
listed in Table fl] 



III. THE PERTURBATIVE QCD CALCULATION 

The PQCD approach is an effective theory to handle hadronic B s decays. Because it 
takes into account the transverse momentum of the valence quarks in the hadrons, one 
will encounter double logarithm divergences when the soft and the collinear momenta 
overlap. Fortunately, these large double logarithm can be re-summed into the Sudakov 
factor [261]. There are also another type of double logarithms which arise from the loop 
corrections to the weak decay vertex. These double logarithms can also be re-summed and 



resulted in the threshold factor [27J . This factor decreases faster than any other power of 
the momentum fraction in the threshold region, which removes the endpoint singularity. 
This factor is often parameterized into a simple form which is independent on channels, 



twists and flavors [28|. Certainly, when the higher order diagrams only suffer from soft 



or collinear infrared divergence, it is ease to cure by using the eikonal approximation 



29J. Controlling these kinds of divergences reasonably makes the PQCD approach more 
self-consistent. 

Here we take the decay B° s — > a\K° as an example, whose all of the single hard gluon 
exchange diagrams are shown in Figure 1. These diagrams contain all of the leading 
order contributions to the decay B° s — > a\K° in the PQCD approach. Diagrams 1(a) and 
1(b) are called factorizable emission diagrams, where ai(1250) is at emitted position, the 
corresponding amplitude is presented as F c k- In the PQCD approach, the form factor can 
be extracted from this amplitude. If ai(1250) is replaced with &!(1235), the amplitude 
contributed by the (V — A)(V ± A) operators would be zero due to the vanishing decay 
constant f^. Diagrams 1(c) and 1(d) are called nonfactorizable emission diagrams, the 
corresponding amplitude is represented as M ek . For the decay B° s — > a®K°, F ek included 
the color suppressed Wilson coefficient C2/S + C1 gives the dominated contributions, while 
for the decay B° s — ¥ b®K , M ek included the color allowed Wilson coefficient C2 gives the 
dominated contribution. Diagrams 1(e), 1(f) and l.(g), l-(h) are called nonfactorizable 
and factorizable annihilation diagrams, respectively, and the corresponding amplitudes 
are written as M ak and F ak . For the decays B° s — )■ a 1 (6 1 )7r, only these annihilation type 
amplitudes can contribute to the final results. If the meson K° is replaced with a vector 
meson K*(p,u,(p), the amplitudes will become complicated, for both longitudinal and 
transverse polarizations can contribute to the decay width. So we can get three kinds 
of polarization amplitudes Ml (longitudinal) and Mn,t (transverse) by calculating these 
diagrams. Because of the aforementioned distribution amplitudes of the axial-vectors 
having the same format as those of the vectors except a factor, so the formulas of here 
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FIG. 1: Diagrams contributing to the decay B® — > a®K°. 

considered £?° — ¥ ai(&i)V decays can be obtained from the ones of B° s — > VV decays by 
some replacements. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 



We use the following input parameters in the numerical calculations [30|, [31 



f Bs = 230MeV, M Bs = 5.37GeV, M w = 80.41GW, 
r Bs = 1.472 x l(r 12 s,a = 91.0 o ,7 = 67.2 o , 
\Vtd\ = 8.58 x 1(T 3 , \Vts\ = 0.03996, \V tb \ = 0.999, 



W, 



ud 



0.97425, IK 



0.22539, \V ub \ = 3.54 x 10' 



(16) 

(17) 
(18) 
(19) 



In the -Barest frame, the decay rates of B s — > ai(6i) V, where V represents K*, p, u, (ft, 
can be written as 



n 2 d -r 2 ) 



32nM B 



(20) 



a=L,N,T 



where M° is the total decay amplitude of each considered decay. The subscript a is 
the helicity states of the two final mesons with one longitudinal component and two 



transverse ones. The decay amplitude can be decomposed into three scalar amplitudes 
a, b, c according to 
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'a)P^R 



2 7 - r 3p? 



(21) 



where M 2 and M 3 are the masses of the two final mesons ai(&i) and K*(p,u,(p), respec- 
tively. The amplitudes A4l, -M-n, -M-t can De expressed as 



M L = a e* 2 (L) ■ e* 3 (L) + 
M N = a, M T 



Ml 



M 2 M 3 



M 2 M 3 
c. 



e* 2 (L) ■ P 3 e* 3 (L) ■ P 2 , 



(22) 



We can use the amplitudes with different Lorentz structures to define the helicity ampli- 
tudes, one longitudinal amplitudes H Q and two transverse amplitudes H±: 



H = MlM L , H± = M 2 B M N =F M 2 M 3 Vr 2 - 1M T , (23) 

where the ratio r = P 2 -P 3 /(M 2 M 3 ). After the helicity summation, we can get the relation 

J2 M aj M a = \M L \ 2 + 2(\M N \ 2 + \M T \ 2 ) = \H \ 2 + \H + \ 2 + \H_\ 2 . (24) 



a=L,N,T 



The matrix elements Aij of the operators in the weak Hamilitonian can be calculated 
by using PQCD approach, which are written as as 



Mj - V ub V* d ( s) Tj - V tb V t * d ^Pj 

= K^^T^l + ^e^)^)), 



(25) 



where j = L,N,T and a and 7 are the Cabibbo-Kobayashi-Maskawa weak phase angles, 
defined via a = arg[— v tb yt ] and 7 = arg[— v *ly* ], respecitvely. Sj is the relative strong 



phase between the tree and the penguin amplitudes, which are denoted as "T/' and "Pj" , 
respectively. The term Zj describes the ratio of penguin to tree contributions and is 
defined as 



V tb V*, 



td(s) 



v ub vi 



ud(s) 



P, 



T 



(26) 



In the same way, it is easy to write decay amplitude Aij for the corresponding conjugated 
decay mode: 



M 



t - V ub V ud(s)Tj - V t lV td ( s) Pj 



V* b V ud{s) Tj(l + Zj e 



,J(-a(-y)+5j)\ 



(27) 



So the CP-averaged branching ratio for each considered decay is defined as 

B=(\M S \ 2 + \M J \ 2 )/2 = \V ub V: d{s) \ 2 [Tl(l + 2z LC osa(~f)cos5 L + zl) 



+2 2_. Tj(l + 2zj cos a (7) cos 5j + z? 



j=N,T 



(28) 



Like the decays B® —¥ VV, there are also 3 types of helicity amplitudes, so corresponding 
to 3 types of Zj and Sj, respectively. Compared with the decays B® -¥ ax(bi)V, the cal- 
culation formula for the branching ratios of other considered decay modes B® — > ai(bi)P 
are simpler, for only the longitudinal polarized component of the axial-vector combining 
with the distribution amplitudes of the pseudo-scalar meson can contribute to the final 
branching ratio. 

Using the input parameters and the wave functions as specified in this section and 
SecJTTl it is easy to get the branching ratios for the considered decays which are listed 
in Table [TTJ where the first two errors are the B s wave function shape parameter u;& = 
0.5 ± 0.05 GeV and the B s meson decay constant fs s = 0.23 ± 0.02 GeV, respectively. 
The third error is induced by the hard scale-dependent varying from A.q CD = 0.25 ± 0.05 
GeV. The last error is from threshold resummation parameter c, varying from 0.3 to 
0.4. The dominant topologies contributing to these decays are also indicated through the 
symbols T(tree), P(penguin), Pg;vK(electroweak penguins), C(color-suppressed tree) and 
ann (annihilation). 



A. B° s -)• ai{bx)K{K*) 

The decays B® — > a^ K + (K* + ) have the contributions from the factorization emission 
diagrams with a large Wilson coefficient C2 + Ci/3 (order of 1), so they have the largest 
branching ratios and arrive at 10 -5 order. While for the decays B° s — > aiK°(K*°), the 
Wilson coefficient is C1 + C2/3 in tree level and color suppressed, so their branching ratios 
are small and fall in the order of 10~ 7 ~ 10 -8 . Although the decay B° s — > a\K° is tree 
dominated, the contributions from tree operators between the factorization and nonfac- 
torization emission diagrams cancel each other mostly, which induces its tree amplitudes 
to have a very small real part. It does not happen in the channel B° s — V a^K* . At 
the same time, there exist three polarization states for the final mesons and the trans- 
verse polarizations are about 30%. So the decay mode a\K*° has a larger branching 
ratio compared with the mode a\K°. For the decay B® — > b^K , the amplitude of the 
nonfactorization emission diagrams Mj k (T denotes the contribution from tree operators) 
including the large Wilson coefficient C2 receives a larger value, which is about 5 times the 
decay a^K . Furthermore, because of the vanishing decay constant f^, the amplitude F e x 
becomes zero for the decay b^K , while which has large value but the opposite sign with 
amplitude M e k for the decay a\K Q . So one can find that there is much larger contribu- 
tion from the tree operator for the decay b\K° than that for the decay a\K° . The decay 
B° s — > b\K* Q has large branching ratio, which is also because of the large contribution 
from the nonfactorizable emission diagrams. 



TABLE II: Branching ratios (in units of 10 -6 ) for the decays B® — > ai(bi)K(n, n, rf) and B® — >• 
ai(bi)K*(p,uj,<p). In our results, the errors for these entries correspond to the uncertainties 
from the B s meson wave function shape parameter cjg, the B s meson decay constant /#, , the 

(5) 

QCD scale A.qq D and the threshold resummation parameter c, respectively. 





Class 


Br(10 -6 ) 




Class 


Br(l(r 6 ) 


B° s -+ a\K° 

B° -> a^K+ 
5° -»■ a^vr+ 
5° -+ a|vr- 
4° -* a°vr 

B° -»■ a??? 

4° -+ a?V 


C 

T 
ann 
ann 
ann 
Pew 
Pew 


n n«i+0-0 16 +0-°o 5 + -0 13 +0- 029 
U.UO±_ .010-0.005-0.011-0.029 

91 a+S. 1+0.1+0. 9+7.0 

zl -^-5. 5-0.0-1. 5-7.0 

9 7+0.7+0.2+0.3+0.0 

z -' -0.5-0.1-0.4-0.0 

1 Q+0.5+0.0+0.2+0.0 

llO -0.4-0. 1-0.3-0.0 

9 9+0.7+0.1+0.4+0.0 

z - z -0.4-0.0-0.2-0.0 

n -1 9+0.04+0.00+0.00+0.03 

u -- LZ -0. 04-0. 00-0. 02-0.03 

fl qn+0-09+0. 02+0. 00+0. 10 

u - ou -0. 08-0. 01-0. 03-0. 10 


B° s -> alK* 

B° s -+ a^K* + 

B° s -> a^p+ 

B° s -> a\p- 

B° s ->• a lP 

B° s -+ a\u 

B° -»• a?0 


C 

T 

ann 
ann 
ann 
ann 
Pew 


fl RQ+0- 19+0.03+0. 10+0. 12 
u - Dy -0. 13-0.04-0. 12-0. 12 
9 q a + 10. 3+0.1+0. 6+9.8 
zy -*-7.2-0. 1-1.8-9.8 
n qo+0.3+0.1+0.5+0.8 
u - oo -0.3-0.1-0.7-0.8 
n q 7 +0. 2+0.1+0.3+0.4 
u - ' -0.5-0.1-0.7-0.4 
n qo+0. 3+0.1+0.5+0.6 
u -°°-0.2-0.1-0.6-0.6 

n nn4Q+ -000 3 + -000 3 + -00 05 + -000 4 

U.UU^y_ 0003-0 .0004-0 .0002-0 .0004 

n qq+0. 13+0.00+0.02+0. 12 
u -°°-0. 08-0. 00-0. 03-0. 12 


4° -»• 6°K° 

B° s -»■ 6^tf+ 

B° -»■ 6^+ 

4° -> 6+vr- 

4° -* 6?vr° 

B° -> b°7/ 

B° s -+ b\rf 


T 

C 
ann 
ann 
ann 
Pew 
Pew 


9 q+0. 5+0.1+0. 4+0.1 
z -°-0.4-0.0-0.3-0.1 
1 q+0. 1+0.0+0. 1+0.0 
-"-•°-0. 2-0.1-0. 3-0.0 

n n7Q+0-0 13 +0-0 01 +0-00 6 +0-000 

U.U 1 y_o.013-0. 000-0.004-0. 000 
n -, 7+0.02+0.00+0.02+0.00 
u - L ' -0.02-0.00-0.01-0.00 

n ns^+0-0 25 +oooo+ -00 2 +0-000 

u - uod -0.017-0.000-0.013-0.000 

n -1 q+0. 05+0. 01+0. 01+0. 01 
u - io -0.01-0.01-0.00-0.01 
n q 9 +0. 09+0. 02+0. 00+0.02 
u - oz -0. 04-0. 00-0. 01-0. 02 


4° -> &$tf*° 
£ s ° -+ 6+p- 

4° -»• &?</> 


T 

C 
ann 
ann 
ann 
ann 
Pew 


q c+0. 6+0.1+0. 6+0.1 

°- d -0.5-0.2-0.6-0.1 

9 n +0.2+0. 1+0.2+0.3 

z - u -0.2-0. 1-0.3-0.3 

fl OQ+0. 06+0. 01+0. 19+0.05 

u - oo -0. 08-0. 02-0. 18-0.05 

1 1+0.1+0.0+0.2+0.0 

i - i -0. 1-0.0-0. 3-0.0 

fl nc+0. 04+0. 01+0. 25+0. 03 

u - yd -0. 06-0. 01-0.24-0. 03 

n m 1+0.001+0.000+0.001+0.002 

U.Ull_ ooi-0 .00-0 .000-0 .002 
n 91+0.04+0.00+0.03+0.00 
u - zl -0. 03-0. 00-0. 04-0. 00 



B. B[ 



ai(bi)n(p,u) 



These channels belong to the annihilation type decays, contributed by the 
W— annihilation and W— exchange diagrams. The decays B® -+ ai(bi)n(p) are sensi- 
tive to the wave functions of the final states. If the final mesons are n and a\, the 
branching ratios can arrive at 10~ 6 order, while for the n and b\ final states, the branch- 
ing ratios become 10~ 7 order even smaller. In a word, B(B°. — > Oi7r) > B(B°. -+ bin). 
The condition is contrary for the decay modes ai(bi)p. The branching ratios of decays 
B® -+ p + a>i , p a>i, P^cli are very near each other. There exists the similar case with the de- 
cays B° s — v p + n~ , p°n°, p~n + , whose branching ratios are predicted as (2.2, 2.3, 2.4) x 10~ 7 
[32|, respectively. We also show the Cabibbo-Kobayashi-Maskawa angle 7 dependence of 
the branching ratios of decays B® — > ai(bi)ir(p) in Fig.2. It is easy to see that the branch- 
ing ratio for the decay with two neutral mesons in the final state lies the between those 
of other two decays in most range of < 7 < 180°. 

As for the other two annihilation type decays B° s — > a°_uj, b°u, whose branching ratios 
are in the order of 10 -8 ~ 10~ 9 . It is easy to see that this kind decay is sensitive to 
the quark structure of the final mesons. Compared with the decays B® -> a5(&?)p°, the 
difference is mainly from the signs of dd component in the mesons u and p°, which in- 
duces different interference effects between the amplitudes from the penguin operators: 
constructive for the decays a5(6?)p°, destructive for the decays a (b°)uj. From our calcula- 
tions, we find that the penguin amplitude for the decay di(bi)p° is about 20.4(48.2) times 
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FIG. 2: The dependence of the branching ratios on the Cabibbo-Kobayashi-Maskawa angle 7. 
In these panels, the solid lines are for the decays B® — > a^b^n , a?(6i) p° , dotted lines for 
B® -)■ a±(bi)ir + ,a±(b^)p + , dashed lines for B® —- d£(bf)ir~,af(bi)p~. 



of that for the decay a°(6?)o;. 



C. B 



ai(6i)r?W 



The main contributions to these four decays are from the electro-weak penguin op- 
erators. Although the contributions from the tree operators have a prominent increase 
for the decays B® — > b\ri^' compared with those for the decays B® —¥ a\r]^ . The former 
are about 5(7) times larger than the later. For the tree operator contributions are the 
Cabibbo-Kobayashi-Maskawa suppressed by a factor 50, so the increased tree operator 
contributions for the decays B® — > b\rf'' bring a slight increase to the branching ratios. 

We also checked the sensitivity to the values on the Gegenbauer moments for all the 

considered decays. If one takes smaller Gegenbauer moments, such as af = 0.05 ± 0.02 

33 1 . 0.10 ± 0.12 [3J], a^' = 0.115 [35[, the branching ratios have a few percent change 



for most of decays B® — y ai(bi)7r(K), more than 10 percent change for only very few 
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channels. So we considered that the uncertainties caused by the Gegenbauer moments 
are small and can be neglected. But it is not the case for the decays B° s — > ai(bi)r]('\ If 
one takes the newer Gegenbauer moments as given in Ref. [351 ]: 

^ = 0.115, < = -0.015, (29) 

The branching ratios will have a prominent change, 

B(B° S -+ aln) = (0.97i°;t!i°;K?^:l) x lO" 7 , (30) 

B(B° 8 -* aW) = (2.1j£K£SK!) x 10- 7 , (31) 

B{B» -> fifr) = (0.211g.^^«g) x ID" 7 , (32) 

B(B» -+ b%>) = (0.75+°;KKr^) x 10~ 7 , (33) 

where the errors come from the B s meson wave function shape parameter u>b = 0.5 ±0.05 
GeV, the B s meson decay constant f Bs = 0.23 ±0.02 GeV, the QCD scale h^ CD = 0.25 ± 
0.05 GeV and threshold resummation parameter c varying from 0.3 to 0.4, respectively. 
Especially for the decays B® — y a\r}^'\ their branching ratios are sensitive to Gegenbauer 
moments and increase to 7 ~ 8 times by using the newer Gegenbauer moments. Certainly, 
the increases of the branching ratios for decays B® — > &°^ are n °t so large. It is need to 
clarify which Gegenbauer moments are more reasonable. 

D. S°^ai(6i)^ 

These two decays are dominated by the electro-weak(EW) penguin operators. Though 
their branching ratios are small, these two decays are interesting to invest the effect from 
the electro- weak penguins, where there might exits new physics [36|. The presence of 
a new physics contribution from EW can enhance the branching ratios of the decays 
Bg —$■ 7r(p)0, which are used to improve the B — > nK "puzzle" [37]] - If here considered 
two decays have such effect, it is deserve more research attention. 

V. POLARIZATION FRACTIONS OF THE DECAYS B° -+ ai(&i)V 

For the decays B° s —¥ a 1 (6 1 )l / , another equivalent set of helicity amplitudes are often 
used, that is 

A = -M 2 B M L , 
An = V2M 2 B M N , 



A± = M 2 M 3 V2(r 2 - l)M T . (34) 

Using this set of helicity amplitudes, we can define three polarization fractions /b,||,± : 

I A |2 

fo ^ ~ \A \* + \A l{ \i + \A ± r (35) 

The formalism of the wave function has great influence to the polarization fractions for 
some decays. In Ref. [38| , the author suggested that taking the asymptotic models for the 
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TABLE III: Longitudinal polarization fraction (/^) and two transverse polarization fractions 
(/||) f±) for the decays B® — > ai{b\)V . In our results, the uncertainties of /l,/ii,/j_ come from 
the B s meson wave function shape parameter ujb, the B s meson decay constant /» , the QCD 

(5) 

scale Aqcd an d threshold resummation parameter c, respectively. 





h(%) 


/||(%) 


/+(%) 


B° s -+ a\K*° 

B° s -+ a^K*+ 

S« -+ a]> + 

5° -+ afp- 

B° s -+ a?p° 

4° -+ a?w 

4° -> o?0 


fiS q+6. 1+2.7+1.5+3.9 
°°- y -6.4-2.8-2.4-3.9 

an r+o. 2+0.3+0. 2+0.1 

y "-°-0.3-0.2-0.3-0.1 
q 7 7+0.1+0.5+0.6+0.9 
y '-' -0.3-0.4-1.2-0.9 

Q7 O + 0.2 + 0.4+0.6+1.0 

y '-"-0. 2-0.3-1. 1-1.0 

Q7 O + 0.2 + 0.3 + 0.6+1.0 

y ' -"-0.1-0.3-1.0-1.0 

00 .+1.0+2.4+3.5+6.1 
"°-^-0. 9-2. 2-2. 5-6.1 

Q A Q+0.0+0.1+0.0+0.2 

y ^'"-0. 0-0.1-0. 0-0.2 


ic 1+3.1+1.3+0.9+2.0 
ld - 1 -3.0-l. 2-0.9-2.0 
a q+0.1+0. 0+0.2+0.1 
^• J -0.1-0. 1-0.0-0.1 
9 9+0.2+0.3+1.1+0.8 
^•^-0.1-0.2-0.6-0. 8 
9 1+0.2+1.1+0.3+1.0 
z - 1 -0.2-0.5-0.2-1.0 
9 1+0.1+0.3+1.0+0.9 
z - i -0.1-0.3-0.6-0.9 

Q O + 0.5+1.4+1.0+2.8 

y -"-0.4-1.3-2.2-2.8 
9 O+0.0+0.0+0.0+0.1 
z -"-0.0-0. 0-0.0-0.1 


ifi n+ 3 - 4 + L5 + L6 + 1 - 9 

lu -"-3.1-l. 4-0.8-1. 9 
A C+O.l+0. 1 + 0.2+0.1 
^• d -0.2-0.1-0.1-0.1 

n 1+0.0+0.0+0.1+0.1 

"• -0.0-0.0-0.0-0.1 

n 1+0.0+0.0+0.0+0.1 
"•-•--o.o-o.o-o.o-o.i 
n 1+0.0+0.0+0.1+0.1 
"•-•--o.i-o.o-o.i-o.i 

n n + 0.3+0.8+1.5+2.3 

"•"-0.4-0.9-1.4-2.3 
9 .+0.0+0.0+0.0+0.1 
z -*-0.0-0. 0-0.1-0.1 


4° -+ 6°^*° 

5° -+ b^K*+ 

B° s -+ 6+p" 

4° -+ 6^/9+ 

5° -+ 6?p° 
4° ~> bfy 


qq 9+0.2+0.2+0.2+0.3 
y "- z -0. 4-0.2-0.4-0.3 
qa 1+0.7+0.6+0.8+1.7 
y ^- 1 -0. 7-0.6-1. 2-1. 7 
Qfi q+0. 3+0.5+0. 9+2. 7 
yo - y -0.3-0.6-2.3-2.7 
Q1 R +0.4+l. 5+3.1+4.5 
yi -"-0. 6-1.4-5. 7-4.5 
QC r n +0. 2+0.8+1. 9+3.8 
yd -"-0. 4-0.9-4.1-3.8 
R q ,.+3.3+3.7+12.7+12.2 
"°-^-2.6-3.8-12.5-12.2 
QQ c+0. 0+0.0+0. 0+0.0 
yy - d -0. 0-0.0-0. 0-0.0 


n q+0.1+0. 1+0.1+0.1 

"• y -0.1-0. 1-0.1-0.1 

9 n+0.3+0.3+0.6+0.8 

z -°-0.3-0.2-0.4-0.8 

9 q+0.2+0.4+1.8+2.0 

z -°-0.2-0.4-1.0-2.0 

1+0.5+1.4+5.5+5.5 

"•-•--0.4-1.2-3.0-5.5 

a 7+0.3+1.1+3.8+3.6 

*•' -0.2-1.2-1.8-3.6 

91 7+1.7+2.4+7.6+8.2 

Zi -' -2.0-2.2-7.5-8.2 

n 9C+0. 01+0. 00+0. 01+0. 00 

"• zd -0. 03-0. 00-0. 03-0. 00 


n q+0.1+0. 2+0.3+0.1 
"• J -o.i-o. 1-0.1-0.1 
q 1+0.3+0.3+0.7+0.9 
°- 1 -0.4-0.4-0.4-0.9 

n 0+0.1+0.1+0.5+0.7 
"•"-0.1-0.1-0.6-0.7 
n q+0.1+0. 0+0.2+0.1 
u -°-o.o-o. 1-0.1-0.1 

n q+0.0+0. 1+0.3+0.3 
"•°-0.0-0. 1-0.1-0.3 
14 o+l-O+l. 5+5.0+4.2 
l*.0_ 1 2 _l .4-5 .1-4.2 
n 9^+0.01+0.00+0.03+0.00 
u - zd -0.01-0.00-0.01-0.00 



K* meson distribution amplitudes instead of its traditional formalism leads to a smaller 
B —¥ K* form factor (Aq ~ 0.3). The smaller form factor responds to the smaller longitu- 
dinal polarization fraction. Another result is that the strengthened penguin annihilation 
and nonfactorizable contribuitons further bring it down. In the decays B® — y a\(bi)K*, 
we also take the asymptotic models for the K* meson wave functions and only find the de- 
cay mode a\K*° with smaller longitudinal polarization fraction about 70%. If we neglect 
penguin annihilation contribution in the decay B° s — y a^K* , and find that the branching 
ratio changes from 6.9 x 10~ 7 to 5.5 x 10 -7 , while the longitudinal polarization receives 
a larger increase and arrives at 93.1%. If we neglect nonfactorizable contribution, both 
the branching ratio and the polarization fractions will become much smaller. Compared 
with B® — y a\K*° and B° s — y b\K*° decays, we argue that the polarization fractions are 
also connected with the symmetric properties of a\ and b\ distribution amplitudes, which 
might have a sensitive effect in the penguin annihilation contribution. If one neglects pen- 
guin annihilation contribution in the decay B® —y b\uj, the longitude fraction can amount 
to 95.4% and the branching ratio decreases by 30%. In a word, the contributions from 
the penguin annihilation diagrams are very sensitive to the final polarization fractions for 
some decays. 

In Table 1111} we list the longitudinal polarization fraction (/x,) and the transverse 
polarization fractions (f\\, f±) for the decays B® —y Oi(bi)V, where the errors come from 
the B s meson wave function shape parameter Ub = 0.5 ± 0.05 GeV, the B s meson decay 
constant f Bs = 0.23±0.02 GeV, the QCD scale A^ D = 0.25±0.05 GeV and the threshold 
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TABLE IV: Direct CP-violating asymmetries (in units of %) for the decays B® 



ai(6i)K(vr,r/,V) and 5° ->• oi(6i)-PC*(p,a;,</>) (except B° -)■ a\K 



t>-*0 „o 



a\uj 



b\uj). In our results, 
i(5) 



the errors for these entries correspond to the uncertainties from w&,/b s , the QCD scale ^qqd 
and the threshold resummation parameter c, respectively. 





Class 


Br(10" 6 ) 




Class 


Br(10 -6 ) 


B° s ->• ar^ + 
£ s ° -»■ c^7T+ 

B s ° -)■ a?vr° 
5° -)• a\n 
B Q S -> a?V 


C 

T 
ann 
ann 
ann 
Pew 
Pew 


fifi -1+10.4+3.3+4.2+37.3 
°°- 1 -6.5-3.3-5.1-37.3 
Q 7+1.4+0.8+0.4+0.7 
y -' -1.6-0.9-0.2-0.7 
9n c + 1. 3+0. 3+0.4+0.2 
zu -°-l. 3-0. 0-0.6-0. 2 
q 9+0.3+0.0+0.6+0.2 
°' z -0. 3-0.1-0. 8-0.2 
-i A n +l. 1+0. 1+0.2+0.0 
J. i i.U_ 1 .2-0.0-0.0-0.0 

qi q + 0.0 + 0.3 + 0.2 + 4.1 

ol -°-2.8-0.2-5.2-4.1 
1n 9+1.4+1.2+2.3+2.1 
iu - z -0.0-l. 3-0.4-2.1 


B° s -> b\K° 
B® -)■ b^+ 

5° -4- 6+7T- 

B° s -»• &?r7 

B° -► 6?!/ 


T 
C 

ann 
ann 
ann 
Pew 
Pew 


a-, /|+5. 3+3. 0+2. 1+0. 3 
* i - ft -5.0-3.1-0.8-0.3 
7A 7+8.1+3.3+0.4+2.5 
''*•' -7.3-2.6-3.6-2.5 
10 7+2.3+0.0+1.4+0.0 
iz -' -3.3-0.1-2.6-0.0 
9/1 c+0.7+0. 0+2. 5+0.1 
z ^- d -3.4-0.0-5.3-0.1 

9 q q + 2.8 + 0.1 + 5.0 + 0.2 

zo -°-1.3-0.1-2.8-0.2 

9 c n+0.0+0. 0+0.5+3. 4 
z,J - u -4.0-2. 8-4.8-3. 4 
99 7+0.0+0.0+0.9+2.4 
zz -' -3.6-2.5-7.2-2.4 


5° -»■ a^iP+ 

B° s ->■ a+p- 
4° -► a?p° 


T 

ann 
ann 
ann 
Pew 


ii 1+1.5+1.0+0.7+1.5 
-'--'-•-'--1.7-0.9-0.5-1.5 
a q+0. 6+0.7+1. 6+1. 5 
^•"J-0. 4-0.5-3.3-1. 5 
fi n +2. 1+1. 1+2. 4+3.1 
"•"-0.8-1.0-3.5-3.1 
a R +l. 3+0.7+1. 6+2.1 
^•"-1.5-0.8-2.9-2.1 
n 9+1.4+0.0+1.4+0.9 
"• z -1.4-0. 0-1. 9-0.9 


5° -)- &°K*° 

B° s ->• 6+p- 
4° -> 6?p° 
4° -»> b?0 


T 
C 

ann 
ann 
ann 
Pew 


9 7+4.2+0.3+5.8+3.2 
z -' -3.7-0.2-5.2-3.2 
n on+7-4+0.3+7.5+3.9 
"•°"-7.4-0.2-6.6-3.9 
qi fi +0.2+0. 0+3.6+0.2 
oi -"-0.2-0.1-2.8-0.2 
n q+0.2+0.2+0.5+0.0 
J - o -0.6-0.3-0.4-0.0 
o q+0.2+0. 1+0.8+0.0 

°--'-o.o-o.2-o.2-o.o 

n oi +0.32+0.00+0.11+0.00 
"• Oi -0. 15-0.00-0. 12-0.00 



resummation parameter c varying from 0.3 to 0.4, respectively. Except the decays B® — > 
a®K*°, a\<jj, b\oj, the longitudinal polarization fractions of other B® —¥ ai(&i)V decays are 
very large and more than 90%. 



VI. DIRECT CP ASYMMETRY 



Now we turn to the evaluations of the CP-violating asymmetries in PQCD approach. 
In view that most of B° s — >■ ai(&i)V decays have small transverse polarization fractions and 
only about few percent. So we can neglect them in our calculations and the expression 
for the direct CP-violating asymmetries of the decays B® — > ai(6i)V (except B® — > 



0^*0 ^o 



a\K 



a\uo 



b\uS) become simple, which can be got by using Eg. ( 125]) and Eq.(J2 



Adir 



\M\ 2 -\M\< 
\M\ 2 +\M\< 



2z] J sma( r y) sin^ 



2zl cos a (7) cos 8l + z 



(36) 



The direct CP-violating asymmetries for the decays B® — > Oi(6i)P have similar expression. 
Using the input parameters and the wave functions as specified in this section and SecJTTl 
one can calculate the PQCD predictions (in units of 10~ 2 ) for the direct CP-violating 
asymmetries of the considered decays, which are listed in Table IIV] where the errors 
induced by the uncertainties of u b = 0.5 ± 0.05 GeV, f B , = 0.23 ± 0.02 GeV, K^ CD = 
0.25 ± 0.05 GeV and the threshold resummation parameter c varying from 0.3 to 0.4, 
respectively. We find the following points: 



13 




-Q 

< 



80 100 120 140 

y (degree) 




80 100 120 140 160 18 

y (degree) 




80 100 120 140 160 180 

y (degree) 



— 35 

~ 30 E- 

Cl : 

sT 25 E" 

■- a. 
■" "20 

< 

15 
10 
5 

-5 
-10 
-15 



| l l l | l l l | l l l | l l l | l l l | l l l | ' i i | ' 







60 80 100 120 140 160 180 

y (degree) 



FIG. 3: The dependence of the direct CP-violating asymmetries on the Cabibbo-Kobayashi- 
Maskawa angle 7. In these panels, the solid lines are for the decays B® — > a5(6i)7T° , a?^?) p° , 
dotted lines for B® — > aj~ (6^~ )tt + , aj~ (b^ ) p + , dashed lines for B® — > af(bf)ir~,af(bf)p~. 



Like the decay B° s — > 7r°K°, whose direct CP-asymmetry is more than 40% predicted 
by several methods [3J 



39l |40j, the decays B° s — y a\(b\)K Q also have large direct 



CP-asymmetries. Unlike the channel B° s — > b~[K + , the decay B° s — > a~[K + has a 
smaller direct CP-asymmetry. It is because that though there are near penguin 
amplitudes in theses two decays, the tree amplitude of the latter is about 3 times 
as large as that of the former, and the sine values of their strong phases are close 
to each other. The direct CP-asymmetries in the decays B° s — > ai(b\)K* are small. 

The direct CP-asymmetries of the decays B® — > &°?r are sensitive to taking different 
Gegenbauer moments for rf'>. If we take the newer Gegenbauer moments given in 
Eq. (l29|) . their direct CP-asymmetries will change not only in magnitudes but also 
in signs. 

The decays B° s — > a>i(bi)p except the channel B® — > b±p + have smaller direct CP- 
violating asymmetries compared with the decays B° s —¥ 01(61)77. The direct CP- 
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violating asymmetry for the decay B° s — > b^p + is very sensitive to the tree operator 
contribution from the nonfactorization annihilation diagrams: if we neglect such 
contribution, its branching ratio can increase 14%, while the direct CP-violating 
asymmetry becomes only 1.3%. In Fig.3, we show the dependence of the direct CP- 
violating asymmetries for the decays B° s — Y ai(bi)iv(p) on the Cabibbo-Kobayashi- 
Maskawa angle 7. 

• There only exist factorization and nonfactorizaiton emission diagrams for the decays 
Bg — > ai(bi)<f). The direct CP-violating asymmetries in these two decays are small, 
because the interactions between tree and penguin contributions are small. From our 
calculations, we find the ratios of penguin to tree amplitudes for decays B° s — > ai<j) 
and B® -> bi<f) are about 0.06 and 0.004, respectively. The strong phases penguin 
and tree amplitudes are only 0.15 and 0.026 rad, respectively. 

• Compared with decays B Q S — > a\{bi)P, most of B® — > ai{b\)V decays have smaller 
direct CP-violating asymmetries. 

VII. CONCLUSION 

In this paper, by using the decay constants and the light-cone distribution amplitudes 
derived from QCD sum-rule method, we research the decays B° s — > ai(bi)P,ai(bi)V in 
PQCD approach and find that 

• The decays B Q S — > a^ K + (K* + ) have the contributions from the factorization emis- 
sion diagrams with a large Wilson coefficient C2 + Ci/3 (order of 1), so they 
have the largest branching ratios and arrive at 10 -5 order. While for the decays 
B® — > a^K (K*°) , the Wilson coefficient is C\ + C2/S in tree level and color sup- 
pressed, so their branching ratios are small and fall in the order of 10~ 7 ~ 10~ 8 . 
For the decays B® — > biK(K*), all of their branching ratios are of order few times 
10~ 6 . 

• For the pure annihilation type decays B® — > ai{b\)p except the decays B® — > a\Ti 
having large branching ratios of order few times 10~ 6 , the most of them have the 



branching ratios of 10 7 order. The branching ratios of the decays B° s — > 
are the smallest and fall in the order of 10~ 8 ~ 10~ 9 . 



"1V"1 



UJ 



The branching ratios and the direct CP-asymmetries of decays B® — > a1(b®)ri^ are 
very sensitive to take different Gegenbauer moments for ?/'). 

Except for the decays B° s — > a^K* , a^u, b^u, the longitudinal polarization fractions 
of other B® — > ai(b\)V decays are very large and more than 90%. 

Compared with decays B° s — >■ ai(bi)P, most of B Q S — >■ a 1 (6 1 )l / decays have smaller 
direct CP-violating asymmetries. 
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